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Abstract 

We explore the Coleman- Weinberg phase transition in regions outside the valid- 
ity of perturbation theory. For this purpose we study a Euclidean field theory with 
two scalar s and discrete symmetry in four dimensions. The phase diagram is es- 
tablished by a numerical solution of a suitable truncation of exact non-perturbative 
flow equations. We find regions in parameter space where the phase transition (in 
dependence on the mass term) is of the second or the first order, separated by a 
triple point. Our quantitative results for the first order phase transition compare 
well to the standard perturbative Coleman- Weinberg calculation of the effective 
potential. 
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During the last two decades, the perturbative Coleman- Weinberg approach |I] to the 
calculation of the effective potential by a loop expansion has been widely used in studying 
field theories that exhibit spontaneous symmetry breaking. Within the standard model 
and for low masses of the Higgs boson this approach predicts a weakly first order phase 
transition in dependence on the mass parameter. This feature has been used to derive a 
lower bound on the mass of the Higgs boson || . A generalization of the Coleman- Weinberg 
calculation for non-vanishing temperature suggests also a first order phase transition in 
dependence on the temperature, at least if the mass of the Higgs Boson is sufficiently 
small 0. The observation that the baryon asymmetry may have been created during 
the electroweak phase transition in the very early universe |4| considerably renewed the 
interest in this subject ||. Indeed, this hypothesis received much encouragement from the 
fact that the loop expansion for the potential predicts a first order phase transition with 
periods outside thermodynamic equilibrium, which is one of the three main ingredients 
for the generation of net baryon number |J. When the experimental bounds for the 
Higgs mass [0] began to diverge from the theoretical expectations ||, two-Higgs-doublet 
extensions of the standard model entered the scene. They provide additional CP violation 
(which is another requirement for baryogenesis) . Also perturbative calculations of the 
effective potential suggest that the phase transition in these models can be sufficiently 
strongly first order and nevertheless remain compatible with the experimental mass bound 
||. However, it is not clear whether these perturbative methods give the right picture 
of the phase transition, since the same methods fail to give a correct description of the 
second order high temperature phase transition |L2] in models with only one scalar. A 
complete description of the phase transitions in two-scalar models is still lacking. 

Another wide range of application of effective potentials concerns inflationary cosmol- 
ogy. In the model proposed originally by Guth |10| as well as in later models [[II] 



an 



effective potential of the Coleman- Weinberg type has been used to allow for a phase of 
exponential expansion of the early universe. These models often do not specify what the 
underlying theory should be, but they nevertheless assume that a general field theoretical 
framework exists which produces this type of potential. 

The high temperature phase transition can be understood in terms of the temperature 
dependent effective potential U(ip, T) for the scalar field p. The crucial feature of a first 
order transition is the existence of two local minima at tp = and p ^ 0, separated 
by a barrier, in a certain temperature range. Similarly, if the zero temperature theory 
exhibits this feature for a certain range of the mass parameter one speaks about a first 
order transition in dependence on the mass parameter. By continuity in T for T = 0, the 
two phenomena are closely related. The renormalization group improved loop expansion 
is believed to work well near the absolute minimum of the effective potential. Already 
in the original paper JT], however, it was pointed out that the description of possible 
additional relative minima is not within the range of applicability of perturbation theory. 
In the phase with spontaneous symmetry breaking the validity of the one-loop results for 
the form of the effective potential at the origin (p = 0) has not been established so far. 
However, a reliable description of both the absolute and the relative minimum is crucial 
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for a firm establishment of a first order phase transition. 

In this letter we perform a non-perturbative analysis which is capable of dealing with 
the above problem. It confirms the picture arising from the results of the Coleman- 
Weinberg loop expansion for the effective potential in the vicinity of the phase transition. 
We employ a new method based on the effective average action [13] which is the analogue 
of the block spin action |14| for continuous space. The dependence of the effective average 
action Tk on the average scale k (one averages over a volume ~ k~ d ) is given by an exact 
evolution equation |L5] ^ In the limit k — > 0, T k becomes the usual effective action 
(the generating functional of the 1PI Greens functions). For k > the effective average 
potential Uk (the non-derivative part of T^) is not necessarily convex. This permits the 
study of several minima of different depth. The effective average potential becomes convex 



in the limit k — > [17| 



We want to solve the flow equation which describes how the effective average potential 
changes its shape as k decreases from some high scale (UV cutoff) to zero. For this 
purpose we first employ a truncation of the exact non-perturbative evolution equation 
by neglecting wave function renormalization effects and ip 6 and higher effective vertices. 
Even in this approximation our approach goes beyond the usual gap equation for the 
mass term since it effectively contains a similar self-coupled equation for the quartic 
scalar coupling. It also properly accounts for mass effects in the scale dependence of 1PI 
vertices and describes how particles with mass m much larger than k decouple from the 
running of the couplings with varying k. These features are not included in the Coleman- 
Weinberg one-loop calculation. The standard perturbative result for the /3-functions of 
the quartic scalar couplings is recovered from our results for small couplings and vanishing 
masses. In a second step we improve the truncation by including the effects of <p 6 and 
if 8 couplings. This permits the investigation of potentials with a rich variety of shapes. 
For example, the potential at large k may have only one minimum (up to the degeneracy 
due to symmetry), whereas at lower k new relative minima may appear as a result of 
integrating out quantum modes with momenta larger than k. Following the evolution of 
the different extrema permits the determination of the nature of the phase transition. The 
reliability of our non-perturbative method for second order phase transitions has recently 
been demonstrated by computing the high temperature phase transition of the O(N)- 
symmetric scalar theory and the corresponding critical exponents to a high precision |12|] . 
Here we address a model with a first order phase transition. 

We consider a simple scalar field theory with two scalar degrees of freedom in four 
dimensions. We limit ourselves to the zero temperature case and study the phase transi- 
tion with varying mass parameter. After discussing the tree level potential we first recall 
the results of the Coleman- Weinberg loop expansion. Then the method of the effective 
average action is briefly introduced and the formalism for the model under investigation 
is given. We study spontaneous symmetry breaking in a model of two scalar fields tpi, <^ 2 



2 See ref. fig] for other version of exact renormalization group equations. 
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with a classical potential symmetric under tp\ — > — fi, y?2 — > — <p2, <Pi — * ^2 

V = -p 2 (p 1 +p 2 ) + ^(p 2 1 +p 2 2 )+g Pl p 2 , (1) 

where pi = |<pf and p 2 = ^V 9 !- This potential is bounded from below provided g > 
—A, A > 0. It is also convenient to use the parametrization 
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V = -p (pi + P2) + A 



2 (pi + P2) 2 + x P1P2 



(2) 



with x = (g — A)/A measuring the deviation from the case with 0(2) symmetry discussed 



earlier with the method of the effective average action [I3[ |12"|| . The structure of this 
potential is described by three different phases. For p 2 < the classical theory is in the 
symmetric phase with the minimum lying at the origin (S phase). If p? > 0, one or both of 
the scalar fields develop a vacuum expectation value. One finds the minimum on the axes 
for x > (AX phase) or between the axes for x < (M phase). For x = the potential 
exhibits an additional 0(2) symmetry. Due to the symmetries of the potential, the minima 
on the axes occur symmetrically at (pi = p 2 /A , p 2 = 0) and (pi = , p 2 = p 2 /A). In 
the M phase, the minimum is always located between the axes at pi = p 2 = p 2 /(g + A). 

The renormalization group improved one-loop potential can be parametrized in terms 
of scale dependent couplings A and g. An important ingredient is the coupling g which 
couples the two scalar fields. In a study of the scale dependence of A this coupling appears 
in the renormalization group equation for A (and vice versa). Two interesting scenarios 
might be induced, which seem to indicate first order phase transitions. First, a term ~ g 2 
appears in the /3-function for A, so that A, instead of running to zero for k as in 
the O(N) symmetric <p 4 -theory, may reach zero during a finite renormalization interval 
and then be driven to negative values. For p? > (AX phase) this would correspond 
to turning the minimum of the effective potential into a saddlepoint. This situation is 
very similar to the scale dependence of the quartic scalar coupling in the Abelian Higgs 
model originally studied by Coleman and Weinberg |IJ (with g playing the role of e 2 ) and 
is generally believed to indicate a first order phase transition. A second interesting case 
occurs for negative couplings g. When g becomes smaller than —A, again the minimum 
flips over. Let us see more quantitatively what standard perturbation theory tells us. The 
one-loop renormalization group equation for the quartic coupling A reads 

' /A ' (9A 2 + /), (3) 



dt 16tt 2 

with a logarithmic renormalization scale t = Ink. Indeed, If g 2 is large enough, A runs to 
zero within a finite interval of t. Since the evolution of the mixed coupling g turns out to 

this interesting feature always shows up for a sufficiently large starting value go/Ao- In 
the case of negative couplings g we are interested in what happens for x approaching —2 
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or g 



-A. Near x 



-2 one finds 



dx 
~dt 



A 
2^ 



(5) 



which, for large enough A and initial value of x close to —2, leads to values of x < —2 
thus flipping the minimum into a saddlepoint! In the same region one has 



dX 
It 



5A2 

8^' 



(6) 



ensuring that there is always a natural parameter range where this phenomenon occurs, 
indicating again a first order phase transition. 

However, some questions remain open. The above renormalization group equations 
do not depend on the masses of the theory and are in this respect "scale invariant" . This 
is a good approximation for mass scales much larger than the characteristic masses of the 
theory. For a first order phase transition the masses do not vanish even for the critical 
parameters. For scales below the critical mass the running of the dimensionless couplings 
should stop. For a precise understanding of the character of the phase transition we want 
to take these "threshold effects" into account. We therefore introduce a scale k which 
acts as an infrared cutoff for the quantum fluctuations, independently of all other possible 
infrared cutoffs such as field dependent masses and external momenta. This replaces the 
effective potential U by an "effective average potential" in which only field modes with 
momenta larger than k are integrated over. 

We start with the evolution equation |L5] describing the dependence of the effective 
average potential on the average scale k(t = Ink) in arbitrary dimensions d 



dUk 
dt 



(pi,P2) =Vd I dx x^ 1 



dP 



+ 



dt [P + Mf P + Mi 



(7) 



Here the variable x denotes momentum squared and v d 1 = 2 d+1 7T2r(|). The inverse 
average propagator 

P 



x 



1-exp (-f) 



contains an effective infrared cutoff for the modes with x < k 2 . The eigenvalues of the 
mass matrix 



Ml 2 = ^ Ul + U *) + U nft + U22 P 2 ± V ~ U2) + Ullpl ~ U22p2 \ + AU u plP2 (9) 

are determined by the partial derivatives of the effective average potential U\ = dUk/dpi, 
U12 = d 2 Uk/dpidp2 etc. We note the close resemblance of the flow equation (0) to a 
renormalization group improved one-loop calculation with infrared cutoff (for k = one 
has P = x). Up to neglecting effects from the wave function renormalization this flow 
equation is an exact non-perturbative equation |15|, 12 1. The appearance of renormalized 



and p-dependent mass terms instead of bare mass terms turns eq. (0) into a self-coupled 
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equation with a structure similar to the gap equation for the renormalized mass. A gap 
equation for the masses can indeed be obtained by computing dMf 2 /dt using partial 
derivatives of eq. (|7|) with respect to p\ and pi- Similar self-coupled equations for the 
quartic couplings follow from taking second partial derivatives of eq. (0). After perform- 
ing the momentum integration the evolution equation (|7|) becomes a partial differential 
equation for £4(pi, P2, t)- This is difficult to solve explicitly and we will first discuss an 
approximate solution based on neglecting partial derivatives of higher than second order, 
i.e. terms ~ Um or (p & couplings and so on. This approximation is too crude to account 
well for complex structures of Uk, as for example several inequivalent local minima (see 
below) . Nevertheless it is a valid approximation as long as k is much larger than the phys- 
ical masses. Before turning to less restrictive truncations we use the quartic potential here 
as an illustration of our method. The evolution of the potential is then described by a 
set of three renormalized parameters which depend on the average scale k. We choose 
A = Un, x = Uu/Un — 1 where partial derivatives are evaluated at the minimum, and 
a third parameter which is either a mass m? = U\ (S phase) or the distance po of the 
minimum from the origin, i.e. p = (pi) min (AX phase) or p = (pi) min + (p 2 ) m in (M 
phase). At the ultraviolet cutoff A, the physics is described by the classical potential V, 
while at some lower scale k < A by the effective average potential Uk which is obtained 
from eq. fl7|) with the initial condition U\ = V. The effective average potential inter- 
polates between the classical potential (k = A) and the effective potential (k = 0) |i5| . 
A solution of the evolution equation for k — > then leads to the effective 1PI vertices 
at zero momentum. The effective average potential does not in general stay a quartic 
potential but receives higher order terms from the renormalization group flow. Near a 
unique minimum, however, it can be taken to a good approximation to be of the same 
form as in eq. fl2|) with renormalized parameters p 2 , A, and x. In this case our truncation 
of higher order couplings (like <p 6 terms) should be a valid approximation. The evolution 
of the renormalized parameters with changing scale k is derived from eq. (0) by partial 
differentiation with respect to p 1 , p 2 Q. 

We next present the evolution equations for the three relevant couplings for the S 
regime, where the partial derivatives defining the running parameters are evaluated at 
the minimum which lies at the origin (p 10 = P20 — 0). This is always the case for 
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-p 2 > 0. One obtains 



^ = -2v d k d -\4 + x)Xlisi^pj (10) 
^ =2v d k d ~\l0 + 2x + x 2 ) A 2 l d 2 s d 2 (jpj (11) 



3 Note that in the AX phase the evolution equations for A = (7u would differ from the one defined as 
A = U22 by higher order terms in powers of A and x. This is due to the choice of a particular truncation 
of the effective average potential Uk which we fix by defining A always through partial derivatives in the 
direction of the minimum as viewed from the origin. 
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^ = - 2v d k d ~\x + l)x(x - 2) A l d s d 2 (jpj , (12) 



where If and l d are constants of order one defined by 



,,/ _ n ,2n-d f j„ i-l^P D -n-l 



£ = y dxx-^—p- n -\ (is) 







The "threshold functions" ( p-j are given by 



and depend on the dimensionless ratio w/k 2 . They approach unity for vanishing argument, 
vanish for large arguments and describe the decoupling of particles with mass greater than 
k. More details can be found in refs. |L3|, [L2||. Another useful equation, though redundant, 
is the running of the coupling g 

% = 2v d k d - 4 (l0 + Ux + Ax 2 ) \ 2 l d 2 s d 2 ( ^) (15) 
at \ k 2 J 

with g = (1 + x)X. We observe that for d = 4 and m 2 = we recover from eqs. (|i~T| ) and 
(p!5|) the standard one- loop /3-functions of eqs. (|]) and (P (making use of t> 4 = 1/327T 2 



and l\ = 1). 

The only non-perturbative content in the approximation of eqs. ( p]) and (p!5|) arises 
through the threshold functions S2 which describe in a natural way that particles with 
m 2 ^> k 2 should not influence the variation of the couplings with an infrared scale k. The 
presence of mass thresholds in the /3-functions, though physically very reasonable, is not 
seen in many versions of the renormalization group equations. In our approach it arises 
naturally from eq. (|?p, together with an additional equation ( |IDD for the scale dependence 
of the mass term. The "quadratic renormalization" of the mass term given by eq. (pTD| ) 
incorporates in the renormalization group framework the physics related to the "quadratic 
divergences". Although known in practice since a long time in the Wilson approach to 
the renormalization group (for example in lattice studies) this equation is often missing 
in the framework of perturbative renormalization group equations. The reason is simply 
that if no mass scale other than m 2 is present in a given formulation of the renormalized 
potential, ratios like m 2 jk 2 cannot be formed and such approaches are necessarily blind 
to mass thresholds. 

In the AX regime (x > 0) we choose the minimum at (pio = po > , P20 = 0). Here po 
is determined by Ui(p , 0) = and one has U 2 (po, 0) = (g — A)p = x\p . The evolution 
equations now refer to the couplings defined at (po, 0) 

d PW _0„. A,d~2 id 



dt 



2v d k d - 2 If 3 sf(2\K) + (1 + x) sf{x\n) (16) 
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— =2v d k d ~ 4 4 A 2 [9 s d 2 (2Xn) + (1 + xf s d 2 (xXn) 



dx 
~dt 



^ =2v d k d ~ 4 t X 



x(l + x) 



1-2 
x 2 



9 s£(2Ak) + ( y - 2x - 7 ) si(xXn) 



dg 
dt 



--2v d k d - 4 l d A 2 



O 1 27 T 1 9 r 2 



si(2\K) + (1 + 2) 



1 - fx - 2x< 



s d (xXn) 



(17) 
(18) 

(19) 



with k = po/k 2 . If finally — 2 < x < and both scalars develop a vacuum expecta- 
tion value (M regime) we use pio = P20 = \po with p determined by £A(|po>|po) = 
fM^Pcb \po) — 0. The mass eigenvalues are given by M 2 = (2 + x)Ap , Mf = —xXp and 
the evolution equations read 



^ =2^ 2 Z? 
at 



3s?((2 + x)Xk) + - — - sfi-xXn) 



dX 



- -2v d k d ~* ll — L_± 
dt k 1 + x 



3xX 1 + 



2 + x 



s(({2 + x)Xk) -sf(-xXK) 



(20) 



+2v d k d - A l d X 2 



9(1 + -J ^((2 + x)Ak) + (1-~) ^(-xAk) 



:2ii 



ax 



at k 1 + x \ 4 



3 2 + x 



ar 



sf({2 + x)Xk) - s({-xXk) 



-2v d k d ~ 4 l d xX 



9(l + |) 2 s d ((2 + x)XK) + (l-^J 2 s d 2 (-xX K ) 



[22] 



dg 



3xA 1 + f 



^ = - 2^" 4 If — 
at k 1 + x 



s?((2 + x)Ak) -s^(-xAk) 



+2^A; d - 4 A 2 



9(1 + |) 2 4(2 + x)A K ) + (l-|) 2 ^(-zAk) 



(23) 



Again the perturbative one-loop equations of A and g are recovered for d = 4 by putting 
s 2 = 1 in eqs. (|T7|) , ( |l~9|) or eqs. (|2~TD , (^) and non-perturbative effects arise only through 
the mass dependence of the threshold functions. 

With the three sets of evolution equations we are now prepared to study the phase 
transition in the two-scalar model in four dimensions. Fig. 1 shows the flow of the coupling 
x and the dimensionless ratio k = po/k 2 which indicates the position of the minimum of 
the potential. Trajectories are shown with decreasing scale t = ln(fc/A), according to the 
sets of evolution equations given above with initial conditions set at some high momentum 
scale A. The phase diagram of fig. 1 is obtained for A(A) = 0.1. We observe the three 
phases S, AX and M separated by (solid) phase transition lines. The phase transition 
between the AX and M phase occurs always for x = 0. This is easy to understand since 
x = corresponds to an enhanced 0(2) symmetry. This symmetry is preserved if we start 
at A with an 0(2) symmetric potential. Since no trajectory can cross the line x = the 
phase diagram actually splits into two separate parts for x > and x < 0. The phase 
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transition between the AX and M phases as a function of x(A) is of the second order and 
we observe that trajectories with small x are attracted towards the transition line. 

The phase transition between the S phase and SSB phase (AX or M phase) does not 
occur for k = as in the tree approximation but rather for positive k(A). In the S phase 
near the phase transition a typical trajectory starts at A near the critical line at k(A) > 0. 
It stays near the critical line for a certain range of scales, moving according to the arrows 
as dictated by the evolution equations in the AX or M phase. At a certain scale k c it 
deviates strongly from the critical line and k reaches zero at some scale k s > 0. Then 
the evolution equation of the symmetric regime has to be used, with boundary condition 
m2 {k s ) — 0. The true phase diagram is, of course, three-dimensional and we only show a 
projection. The trajectories also have a component perpendicular to the projected plane. 
Since for small A the evolution of A is very slow, the projection on the plane with constant 
A gives a satisfactory picture. 

On the critical line one observes three fixpoints: The "Ising fixpoint" at x = — 1 
corresponds to two disconnected scalar theories. The "Heisenberg fixpoint" at x = is 
characterized by an additional 0(2) symmetry, and the "cubic fixpoint" is situated near 
x = 2 (compare with eq. (|T2"D) [JJ The three points are partial fixpoints and become full 
fixpoints only in the limit where the running of A is neglected. (The only true full fixpoint 
is the Gaussian fixpoint for A = 0.) They become genuine fixpoints in three (or 4 — e) 
dimensions and we have used here their names common in statistical mechanics |18| . The 
fixpoints at x — — 1 and x = 2 are repulsive and the 0(2) symmetric one at x = is 
attractive. 

In the region —1 < x < 2 the phase transition is of the second order. Near-critical 
trajectories stay very near the phase transition line before k deviates either towards zero 
(S phase) or infinity (SSB phase) at some scale k c . Exactly on the phase transition the 
trajectories stay on the critical line for infinitely long "time", i.e. k c — > 0. For x < — 1 the 
critical trajectories run towards x = —2 where the quartic polynomial is unstable. We 
will show that in this region the phase transition is of the first order. The point x = — 1 
is a tricritical point separating the first and second order part of the phase diagram. 
Similarly, for x > 2 the trajectories run to x — > oo (A — > 0). This part of the phase 
diagram corresponds again to a first order transition with tricritical point at x = 2. For 
the establishment of the first order character of the phase transition the two interesting 
regions are x ~ — 2 and x — > oo. They lie on the left and right edges of the phase diagram 
in fig. 1. 

In the regions x ~ —2 and x — > oo (A ~ 0) the approximation of the effective average 
potential £4(pi, p 2 ) by a quadratic polynomial in pi and p 2 becomes inadequate. Indeed, 
if one would insist on solving the corresponding evolution equations discussed before, 
one would find in these regions an asymptotic behaviour of the phase transition line 
characterized by fixpoints in n 3 \(g + A) and K 3 Xg respectively. They simulate a second 
order transition and are artefacts of an insufficient truncation. One has to account for 

4 The true fixpoint occurs exactly for x = 2 and the small deviation from this value is due to a 
particularity of our truncation. For details see ref . J2(]] . 
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the possibility that the effective average potential has more than one local minima in 
these regions. In this case a description of the potential in terms of "local variables" 
at the absolute minimum, i.e. a polynomial expansion around this minimum, may be 
insufficient. A more global approach, which permits the simultaneous existence of more 
than one local minima, and monitors their properties as well as the differences in height 
between them, seems more appropriate. 

In the remaining part we will extend the study by including higher derivative couplings 
as well as by monitoring the extrema of the potential simultaneously, in order to to get 
an accurate picture of the phase transition. We restrict ourselves to the region x —> oo 
(A ~ 0) 0. We parametrize the effective average potential Uk by eight parameters 

U k = -p 2 (pi + ft) + g V/ 9 ? + Pi) + 9oPip2 
+ va{p\(>2 + Pipl) + vb{pI + p\) 

+ IapIpI + 1b(p\p2 + PxP%) + lc{p\ + pt)- (24) 

A convenient set of running parameters in this region is 

m 2 Q = ^(0) = U 2 (0) = -p 2 A = U n {0) = U 22 (0) g = U 12 (0) 

p from L r i(p ) = m 2 = U 2 (p ) 

\i = U n (po) A 2 = f/ 22 (p ) g = U 12 {p ). (25) 

Here we assume a vacuum expectation value at (p 10 = po, P20 — 0) and the partial 
derivatives of the potential are evaluated at p\ = and p% = po- The first three parameters 
describe the potential at the origin while the last five specify the minimum at p\ = p . 
The set of evolution equations is now given by 

dm 2 , , , d\n , , don , . 

dpo 6(Po) dm - 2 c t \ , tt t \ d P° 

ir = -jam ^r = ^ o) + f/i2(po) ^ 

-tt =6i(Po) + ^m(po) -JT = ^2(po) + U 122 (p ) 

dt dt dt dt 

f =e«(po) + EW^^, (26) 

where ^ = dUk/dt(pi, p 2 ) is defined in eq. (7), with subscripts again denoting partial 
derivation in the p-directions. The quantities on r.h.s. of the flow equations (^) are 
computed by neglecting terms with more than four derivatives with respect to pi, p 2 . 
The remaining partial derivatives (as for example Uu 22 ) are expressed in terms of the 
parameters defined in eqs. ( f23| ) by use of the polynomial approximation of eq. (P^p. This 
leads to a system of eight coupled non-linear differential equations 0. 

5 The region x ~ — 2 can be mapped onto this region by a field transformation j2(J. 
6 The reduced system of three equations (|l^) , (|l7|) , (|l9|) is obtained by the truncation va = vb = J A — 
1b = = and selecting the couplings po, Xi and g. 
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We have solved this system numerically, starting with a classical potential given by 
eq. (D and /i 2 (A)/A 2 = 2.22 x KT 4 , A(A) = 1/90, g{A) = 1/10. As long as X^k) 
remains larger than about 10~ 4 we see small quantitative but no qualitative differences as 
compared to the reduced system of eqs. (|I1^),(|lTl),(|i9l). The trajectories look very similar 



to the ones of a second order phase transition. Only if we start with k(A) so close to the 
critical value that \i(k) runs to values much smaller than 10 -4 the first order character 
of the transition is revealed. We have depicted the phase diagram for small values of 
Ai in fig. 2, where we show the trajectories in the (Ai,«) plane. The third relevant [] 
coupling g remains almost constant in this part of the phase diagram. When starting the 
trajectories in the lower right corner of fig. 2 all the remaining "irrelevant" parameters 
(777.Q, Ao, go, m 2 , A2) have to be set with very high precision equal to the values they take 
after the running along the (almost) critical trajectory for Xi(k) > 10 -4 . Since this first 
part of the running extends over many orders of magnitude in the scale k and the critical 
trajectory is unstable in the relevant K-direction, care has to be taken to keep track of 
the critical trajectory accurately 0. 

The right and lower corner of the phase diagram of fig. 2 corresponds to an effective 
average potential Uk with only one minimum. On the other hand, if a trajectory crosses 
the dashed line a second minimum develops and the upper left corner of fig. 2 corresponds 
to a potential with two local minima at p = and p = po 7^ 0. To the left of the dashed- 
dotted line the minimum at the origin is deeper than the minimum at po ^ 0. The curve B 
(which asymptotically approaches this line) corresponds to the phase transition line which 
separates the symmetric phase (all trajectories below B) from the phase with spontaneous 
symmetry breaking (all trajectories above B). The SSB phase can be subdivided into 
region I with only one minimum of Uk and region II with minima at p = and p = 
Po 7^ 0, the second one being deeper. The boundary trajectory A separates these regions. 
Similarly, in region III the potential has two local minima with the deeper one at zero. 
Region III is bounded by the spinodal curve C. This is the trajectory where Ai vanishes 
asymptotically for k — > 0. For trajectories below C and above D the minimum at po ^ 
disappears for some scale &2 > 0, where X\(k2) = 0. Finally, for trajectories below the 
boundary D the location of the minimum n(k) reaches zero for some scale k s > with 
Xi(k 2 ) > 0- The behaviour of all trajectories above A and below D resembles closely a 
second order phase transition. The typical behaviour of the first order transition is only 
visible in a narrow range around the critical line in figs. 1 and 2 which corresponds to the 
region between the curves A and D. We should mention at this point that for trajectories 
between A and D we have stopped the running of the couplings at some scale k cnv ^> 
instead of computing the effective potential for k — > 0. The scale k cnv was chosen large 

7 We do not distinguish here between relevant and marginal couplings in the language of statistical 
mechanics. 

8 The inherent numerical uncertainties do not permit to follow the almost critical trajectories in 
one sweep from k — A to k = 0. Algorithms with a "returning" of n at intermediate k had to be 
developed. For the "starting point" for fig. 2 we used values for the parameters approximately given by 
milk 2 = -1.36 x 1(T 4 , A = 6.74 x HT 5 , g a = 4.18 x 1(T 2 , K = p /k 2 = 2.09, m 2 /k 2 = 8.73 x 1(T 2 , Ai = 
6.39 x HT 5 , A 2 = 1.88 x 1CT 5 , g = 4.17 x 1(T 2 . 
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enough so that k 2 +M 2 > \k 2 , and fc 2 +Mf > \k 2 for all pi, p 2 (cf. eqs. (0), (g)), and small 
enough so that the running couplings approached almost constant values. The potential 
with several minima therefore corresponds to the effective average potential Uk cnv , rather 
than the effective potential Uq. In other words, the quantum fluctuations with momenta 
q 2 < k 2 cnv have not been included. The properties of the phase transition can actually 
be understood in terms of the shape of U kcnv - Including the quantum fluctuations with 
q 2 < k 2 nv will lead to a convex effective potential Uq []17| . The physics related to the 
"approach to convexity" is not relevant for our purpose and rather obscures the simple 
picture presented above. 

For k > k cnv we find that the contribution of the ipi modes to the running of 
Un{pi, Pi = 0,k) is very small as compared to the contribution of the tp 2 modes, once 
Xi(k) becomes smaller than 10 -4 . The <p\ modes may therefore be neglected in this region 
0, corresponding to an omission of the first term on the r.h.s. of the flow equation (|7|). 
With i ( 

M 2 ( Pl , p 2 = 0) = U 2 { Pl ) =m 2 + g oPl + -v A p\ + - lB p\ (27) 

and neglecting the /c-dependence of Wq, go,VA,lB we can easily integrate eq. (^). The 
solution for corresponds to a one-loop formula 

Un(pi, 0) =A + v B p x + Tflcpl 



L 2 { 9o+uApi+ h Bpi ) r dxx {(x+u 2 ( P1 ) 
i r°° \ i 

327T 2 Jo [x + U 2 {pi) P kl {x) + U 2 {pi, 



1 



32tt 2 V" u w '2 w Vio \(x + U 2 ( Pl )) 2 (P kl (x) + U 2 ( Pl )) 2 

1 1 



(28) 

Here k x may be chosen so that \\{k\) = 10~ 4 , with P kl = x 1 — exp (— ^ . All 
couplings on the r.h.s. are taken at the scale ki, with u^, u B , ^ B and 7c re-expressed in 
terms of the couplings defined in eqs. (|25|) . Omitting all couplings except Mq, Ao and go, 
we recognize the usual logarithmic behaviour of the Coleman- Weinberg potential 



2 



/':,(/;,.()) ^ An ^ 



32tt 2 



In 



'k\ +mg + g pi\ k\ 



mi + 9oPi J k( + m^ + g p 1 



(29) 



Eq. (p8|) is in agreement with our numerical solution of the evolution equations and 
constitutes a quantitative improvement since no polynomial truncation of Uk{pi,0) is 
imposed for k < ki. 

In conclusion, a combination of the flow equations for X\(k) > 10~ 4 with the "one-loop 
formula" ( p8|) for Xi(k) < 10~ 4 permits a reliable and quantitatively very precise under- 
standing of the first order phase transition. Our non-perturbative method clearly answers 

9 Even though suspected by the smallness of Ai a similar statement cannot be proven within the 
standard perturbation theory since the saddle point appoximation breaks down for small p\. 

10 Since only the fluctuations of (p% are related to the approach to convexity we can replace here k cnv 
with zero. 
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the shortcomings of perturbation theory and firmly establishes the first order character 
of the transition for A (A) < |g(A). The situation is analogous to the abelian Higgs model 
for which a similar non-perturbative investigation has been performed recently |19 |. Our 



method can be easily adapted to non-vanishing temperature or to the three-dimensional 
model relevant for statistical mechanics. In fact, we have presented our equations already 
for arbitrary dimensions d and we emphasize that the usual perturbative infrared diver- 
gences in three dimensions are absent, due to the presence of an infrared cutoff ~ k. For 
a study of the temperature dependence of the effective potential it is sufficient to replace 
the threshold functions of eq. ([14]) by temperature dependent threshold functions fT2|l , 
the modification arising from the discretization of the zero-component of the momentum. 
Combining the method presented here with the study of temperature dependence in ref. 



12] constitutes an excellent basis for a quantitative understanding of the high temperature 



behaviour of the two-scalar model. This study is presented in detail in ref. [20 
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Figures 

Fig. 1 Phase diagram for A(A) = 0.1. 

Fig. 2 First order phase diagram in the region of small Ai. 
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